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PREFACE 


To  evaluate  the  combination  of  third-order  elastic  constants 
of  cubic  crystals,  one  needs  very  accurate  values  of  the  combination 
of  second-order  elastic  constants  K2.  Often  one  can  obtain  numerical 
values  of  second-order  elastic  constants  for  the  sample  of  interest 
from  existing  data.  Due  to  the  error  propagation  and  the  experimental 
uncertainties  at  the  time  of  the  original  measurement,  however,  there 
may  be  uncertainties  in  the  magnitude  of  K2>  An  alternative  way  to 
arrive  at  K2's  would  be  to  measure  them  at  the  time  the  data  are 
taken.  This  measurement  is  quite  straightforward  since  K2's  are 
directly  related  to  the  sound  velocities,  which  can  be  measured.  The 
extra  time  and  effort  spent  in  doing  the  latter  would  be  justified  by 
an  increase  in  accuracy.  The  purpose  of  this  thesis  is  to  study 
these  alternatives  and  compare  error  propagation  in  an  effort  to  arrive 
at  the  most  accurate  way  to  evaluate  1^.  The  analysis  in  this  thesis 
shows  that  the  question  originally  posed  does  not  have  a  unique  answer 
for  all  samples  under  all  conditions.  One  cannot  decide  a  priori 
whether  the  reference  values  of  C.  ■  should  be  used  or  whether  one  should 

*  J 

measure  the  C.j  each  time  one  measures  Cijk'  The  analysis  given,  how¬ 
ever,  suggests  that  one  should  use  reference  data  rather  than  measuring 
each  sample,  if  the  reference  data  are  as  accurate  as  those  of 
McSkimin.  If  such  accurate  data  are  not  available,  one  has  no  choice. 
One  must  measure  the  K2  directly. 
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CHAPTER  I 


INTRODUCTION:  GENERALIZED  DEFINITION  OF 
ELASTIC  CONSTANTS 

The  application  of  a  time  varying  force  on  a  solid  material 
causes  a  deformation  of  the  solid  and  gives  rise  to  stress  waves.  To 
explain  the  behavior  of  a  material  upon  application  of  stress,  several 
general  relations  between  stress  applied  and  the  resulting  strai  i  have 
been  studied.  One  such  relation  was  put  forward  by  Hooke,  according  to 
which,  for  a  linear  elastic  medium  the  stress  applied  to  a  material  is 
proportional  to  the  resulting  strain,  where  the  coefficient  of 
proportionality  is  a  constant  independent  of  stress,  strain  and  their 
time  derivatives.  Although  Hooke's  law  generalized  to  anisotropic 
media  is  quite  appropriate  for  description  of  many  phenomena,  it  fails 
totally  to  explain  nonlinear  phenomena.  Thus,  a  more  general  approach 
is  required.  Such  a  generalized  approach  was  made  by  Murnaghan  (1951) 
who  started  from  the  definition  of  the  energy  of  a  small  volume  sub¬ 
jected  to  a  homogeneous  strain.  This  approach  later  was  applied  to 
crystals  of  cubic  symmetry  by  Holt  and  Ford  (1967). 

If  $  is  the  internal  energy  of  a  unit  mass  of  material  in  an 
undeformed  state,  i.e.,  energy  of  the  solid  in  equilibrium,  then  for 
small  deformations  we  can  expand  the  strain  energy  in  a  series 

A  —  A4-'*+'++/f.  +  (1*1) 

*  0  rl  v2  *3  ’  ' 

where  ^  is  the  first-order  perturbation  term,  etc.  In  terms  of 
elastic  mode1!-  one  can  write  the  same  expansion: 


(1.2) 


>0  +  V.  Cij  nij  +  2!  Ci jk£  nij  nk£ 


+  3!  ^ijk£mn  nij  nk£  nmn  +  •  •  • 


where  there  is  a  summation  over  repeated  indices  which  take  successive 
values  of  1,  2,  and  3,  and  where  C .  - ' s  are  the  elastic  constants  and 

’  s) 

are  the  components  of  a  strain  tensor.  The  first  term  on  the 
right-hand  side  of  Eq.  (1.2),  <t>g,  is  independent  of  strain  and  there¬ 
fore  can  be  set  equal  to  zero  without  loss  of  generality.  The  second 
term,  <j>,,  is  also  set  equal  to  zero,  as  it  corresponds  to  displacement 
without  deformation.  Equation  (1.2)  then  reduces  to: 


2!  Cijk£  nij  n 


k£ 


+  3!  Cijk£mn  nij  nk£  ^mn  + 


(1.3) 


where  are  the  second-order  elastic  constants  and  are  the 

third-order  elastic  constants. 

The  expansion  of  strain  energy  in  terms  of  strains  can  be 
substituted  into  Lagrange's  equation  to  obtain  a  completely  general 
wave  equation  capable  of  describing  both  linear  and  nonlinear  wave 
phenemona  in  solids  of  any  crystalline  symmetry.  Such  an  equation  has 
been  derived  and  has  been  specialized  to  cubic  symmetry.  For  longitudinal 
waves  along  the  principal  directions  in  a  cubic  crystal  the  equation 
takes  the  form 


K  3  “  +'3K  +  * 

6  da6  J 


2 

3__u 

?a2  ja 


(1.4) 


'„*■ 
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where  the  symbol  ^  stands  for  linear  combinations  of  second-order 
elastic  constants  and  K3  stands  for  linear  combinations  of  third-order 
elastic  constants  as  shown  in  Table  1.1.  The  solution  to  the  nonlinear 
wave  equation  is  as  follows: 


'  3K  +  K  1  2  2 

u  =  A i  sin(ka  -  wt)  -  — - -  k  a  cos  2(ka  -  ut) 

^  j 

+  .  .  . 


(1.5) 


This  solution  shows  that  in  a  nonlinear  solid  the  propagation  of  an 

initially  sinusoidal  wave  generates  a  second  harmonic  whose  amplitude 

is  a  linear  function  of  propagation  distance  a  and  is  proportional  to 

3K2  +  K- 

the  combination  of  elastic  constants  - p - .  This  combination  often 

K2 

is  called  the  nonlinearity  parameter.  Measurement  of  the  amplitude  of 
the  second  harmonic  generated  as  a  longitudinal  ultrasonic  wave  propa¬ 
gates  along  the  symmetry  directions  in  cubic  crystals  has  led  to  values 
of  the  third-order  elastic  constants  of  copper,  germanium,  silicon, 
KZnF^,  SrTiO^,  and  CsCdF^,  between  room  temperature  and  77  °K,  or  even 
3  °K. 

But  in  the  course  of  evaluating  the  third-order  elastic  constant, 
K3>  the  experimenter  must  obtain  values  of  the  second-order  elastic 
constants  Kz>  and  the  accuracy  of  the  K3  is  directly  dependent  upon  the 
accuracy  of  K^. 

Since  the  second-order  elastic  constants  of  a  number  of  crystals 
have  been  measured,  it  often  is  possible  to  obtain  numerical  values 
for  the  sample  of  interest  from  existing  data.  But  the  uncertainties 
in  the  magnitude  of  «2  depend  upon  error  propagation  as  well  as  upon 
the  experimental  uncertainties  at  the  time  of  the  original  measurement. 


An  alternative  way  to  arrive  at  expressions  for  Kg  would  be  to 
measure  them  at  the  time  the  data  are  taken.  Since 


this  measurement  is  quite  direct— a  measurement  of  the  density  of  the 
material  and  of  the  velocity  v  of  a  longitudinal  wave  along  the 
direction  for  which  Kg  is  defined.  An  advantage  of  this  procedure  is 
that  Kg  data  are  taken  on  the  same  sample  as  the  data.  But  one 
disadvantage  is  that  such  a  measurement  would  essentially  double  the 
time  spent  in  measuring  the  third-order  elastic  constants  as  a  function 
of  temperature.  Since  such  measurements  typically  require  an 
uninterrupted  period  of  at  least  36  hours,  such  an  increase  in  time 
would  be  justified  only  by  an  obvious  increase  in  accuracy  of  the  data. 

The  purpose  of  the  present  thesis  is  to  evaluate  the  alternatives 
by  comparing  error  propagation.  To  do  this,  measurements  were  made  of 
the  velocity  of  propagation  in  the  three  principal  directions  In 
germanium,  a  cubic  crystal. 

The  data  are  subjected  to  a  detailed  error  analysis  and  the 
results  are  compared  with  the  errors  propagated  when  the  data  of 
McSkimin  (1963)  are  used  to  calculate  velocities  of  longitudinal  ultra¬ 
sonic  waves  in  germanium.  The  comparison  is  made  more  significant  by 
the  fact  that  the  sample  actually  used  by  McSkimin  was  available  for 
the  measurements  reported  in  this  thesis. 


CHAPTER  II 


THEORY  OF  ELASTIC  WAVE  PROPAGATION 

The  evaluation  of  the  second-order  elastic  constants  from  data 
on  the  velocity  of  an  ultrasonic  wave  in  a  cubic  solid  requires  the 
derivation  of  the  wave  equation  from  the  basic  definitions.  The  basic 
definition  of  strain  energy  (Eq.  1.3)  is  general  enough  to  allow  a 
description  of  nonlinear  phenomena;  however,  for  purposes  of  deriving 
the  linear  form  of  the  wave  equation  to  be  used  throughout  the  remainder 
of  this  thesis,  the  strain  energy  can  be  approximated  by  including  only 
the  first  set  of  terms.  Since  these  terms  contain  second  powers  of  the 
strain,  the  coefficients  are  known  as  second-order  elastic  constants. 
Truncation  of  the  strain  energy  expansion  in  this  way  also  allows  one 
to  make  a  definition  of  Hooke's  law  generalized  to  anisotropic  media. 
Although  the  generalized  Hooke's  law  has  limited  value  since  its  very 
definition  prohibits  the  consideration  of  third-order  elastic  constants, 
it  still  is  used  in  engineering  applications  often  enough  to  justify 
its  consideration  in  this  thesis. 

A.  RELATION  BETWEEN  STRESS  AND  STRAIN  IN  AN  ANISOTROPIC 
SOLID  (GENERALIZED  HOOKE'S  LAW) 

To  derive  Hooke's  law  generalized  to  anisotropic  media  It  Is 
adequate  to  retain  only  the  second  power  terms  in  the  strain  energy. 
Equation  (1.3)  thus  becomes 
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C1jk£  nij 


k£ 


) 

(2.1)  I 


Now,  consider  the  work  done  when  a  strain  results  from  a  stress 
o^..  The  work  E  done  by  the  stress  in  causing  the  strain  Is 


E  ’  2  °1J  n1j 


(2.2) 


Equating  the  work  done  to  the  strain  energy  results  in 


°1j  =  C1jk£  \t 


(2.3) 


This  is  known  as  the  generalized  Hooke's  law  for  an  anisotropic  medium. 
In  this  form  it  can  be  used  to  describe  a  medium  of  any  crystalline 
symmetry  in  the  linear  approximation. 


B.  SECOND-ORDER  ELASTIC  CONSTANTS  OF  CUBIC  CRYSTALS 


From  Eq.  (2.3)  It  appears  that  there  must  be  34  =  81  elastic 
constants  C^^.  The  total  number  of  independent  constants,  however, 

is  reduced  by  lattice  symmetries.  When  the  stress  tensor  and  the 
strain  tensor  are  symmetric  (i.e.,  =  o ^  and  nkj,  =  n£k),  the 

number  of  independent  second-order  constants  reduces  to  6  =  36.  This 
reduction  is  done  easily  by  relndexing  the  stress  and  strain  components 
according  to  Voigt  rotation  as  follows: 


Q 

II 

°n 

'  °4  =  °23  =  °32 

°2  = 

°22 

’  °5  =  °31  =  °1 3 

°3  = 

°33 

‘  °6  =  °12  =  °21 

and  (2.8) 


nl  = 

nll 

9 

n4  = 

n23  =  n32 

n2  = 

n22 

9 

n5  = 

n31  =  nl 3 

n3  = 

n33 

9 

n6  = 

n12  =  n21 

With  this  notation  the  linear  relation  between  stress  and  strain  Is 
written  as 

aa  =  CaB  nB  (°’B  =  1§  2>  3>  4>  5»  6>  (2.9) 

where  Cag  are  constants  and  the  range  of  summation  now  is  up  to  6  for 
the  indices  a  and  B.  The  elastic  constants  C  „  and  C-..„  are  related. 

cip  IJKJc 

Some  examples  of  how  these  are  related  are  shown  below: 

C11  =  Cllll 

C1 2  =  C1 122  (2.10) 
C14  =  2C1 1 23  =  2C1 1 32  . 

We  observe  that  Ci^kr  =  Cjikr  and  This  implies  that 

there  are  only  six  combinations  for  the  first  pair  of  indices  and  six 
of  the  second  pair.  Thus,  for  a  cubic  symmetry  since, 


C11  =  C22  =  C33 
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C12  =  C21  “  C13  =  C31  =  C23  "  C32 


C44  =  C55  "  C66 


the  total  number  of  Independent  elastic  constants  reduces  to  three, 


viz.  C 1 1 ,  C -j 2  and  C^. 


C.  RELATION  BETWEEN  C^.'S  AND  SOUND  VELOCITY 


1.  Wave  Speeds  in  an  Elastic  Medium 


The  linear  equation  of  motion  for  a  wave  through  an  elastic  medium 


can  be  derived  by  using  the  strain  energy  Eq.  (2.1)  in  Lagrange’s  equation. 


As  given  by  Green  (1973)  it  is: 


3^uk 

pUi  =  Ci jkJt  dx.^  3Xj 


(2.11) 


where  C..^  are  the  second-order  elastic  constants,  u^  is  the  displace¬ 


ment  at  a  time  t,  and  p  is  the  mass  density  of  the  homogeneous  medium. 


In  the  principal  directions  [100],  [110],  and  [111]  in  a  cubic  lattice 
the  three  equations  represented  In  Eq.  (2.11)  are  uncoupled,  and  to 
derive  a  relation  between  wave  speeds  and  C^’s  in  these  directions  we 
assume  a  solution  to  Eq.  (2.11)  of  the  form: 


u^.t)  -  A0a,  exp  1(ut  -  kmxm) 


(2.12) 


where  AQ  is  the  amplitude  of  the  wave,  ai  are  the  direction  cosines 


of  the  displacement  vector,  km  is  the  wave  vector  given  by 


-  k'm  ’  <T>  • 


(2.13) 


4  *.  \  ' 

'  V.V.TVV, 


m 


where  k  is  the  wave  number,  i  (=  t,m,n)  are  the  direction  cosines  of 

m 

the  normal  to  the  plane  wave,  and  x  is  the  wavelength.  Substituting 
Eq.  (2.12)  into  Eq.  (2.11)  and  using  Eq.  (2.13)  gives  a  set  of  equations 
of  the  form 


tCijk£  ‘  £«.  *j  "  pv  6ik^uk  =  0 


where  v,  the  speed  of  propagation  of  the  wave,  is  given  by 


2  2..  2 

v  =  <jj  /k 


(2.14) 


(2.15) 


and  ui  =  ikuk'  For  a  nontrivia1  solution  of  Eq.  (2.14)  to  exist,  the 
determinant  of  the  coefficient  must  vanish.  Thus, 


I'fjfce  lt  lj  '  ov  sikl  *  0 


(2.16 


For  simplicity,  let  ££  £j  =  X^;  then  (2.16)  reduces  to 


>ik  -  pv  «1kl  =  0  . 


(2.i7: 


Equation  (2.17)  is  an  eigenvalue  problem  of  the  form 

A11  A12  A13  10  0 
A21  A22  A23  -  pv2  0  1  0=0 


A31  A32  A33 


0  0  1 


(2.18; 


whose  solution  gives  the  wave  speeds.  This  determinant  can  be  written 
in  the  form 


“V**  *  -  V  V  V  V  V  V.r/  '/  /  V/Vv 
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x33  ‘  DV 


=  0 


(2.19) 


The  fact  that  this  is  a  3  x  3  determinant  indicates  that  there  are  three 
solutions;  i.e.,  there  are  three  independent  plane  waves,  each  having 
polarization  in  one  of  the  three  orthogonal  directions.  The  speeds  of 
these  waves  are  the  eigenvalues.  To  determine  the  eigenvectors  we  define 
the  direction  cosines  a,  B,  y  of  the  particle  displacements.  They  are 
obtained  by  writing  Eq.  (2.16)  in  the  form; 


11 

X12 

X1 3 

21 

X22 

X23 

31 

X32 

X33. 

have 

used 

5ik  = 

(  > 
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=  PV 


0 
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(2.20) 
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0  1  0 
0  0  1. 

o  repeat,  the  solution  of  the  eigenvalue 


The  solution  of  this  equation 

10  0  1J 

gives  the  eigenvectors, 
problem  (2.17)  for  eigenvalues  and  eigenvectors  gives  information  about 
the  magnitude  and  direction  of  the  speed  of  propagation  of  the  three 
possible  plane  waves  through  a  crystal.  It  is  specialized  to  the  pure 
mode  directions  in  cubic  crystals  in  the  next  section. 


2.  Wave  Speeds  Along  Pure  Mode  Directions  in  Cubic  Crystals 

The  three  principal  directions  in  cubic  crystals  are  shown  in 
Figure  2.1.  In  each  of  these  directions  three  pure  mode  elastic  waves 
(one  longitudinal  and  two  transverse)  may  propagate,  while  in  all  the 
other  directions  waves  may  be  coupled.  We  now  study  the  behavior  of 
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the  elastic  waves  along  the  pure  mode  directions  of  cubic  crystals 


and  derive  relations  between  sound  velocities  and  C^'s. 
a.  [100]  Direction 

For  the  [100]  plane-the-ddrection  cosines  of  the  plane  wave 
normal  are  i  =  1 ,  m  =  0,  n  =  0.  For  this  direction, 

X11  =C11 

X12  =  X13  =  X23  =  X21~X31  =  X32  =  0 
X22  =  X33  =  C44  ' 

Substitution  of  these  values  into  (2.17)  and  simplification  leads  to: 
(Cn  -  pv2)(C44  -  pv2)(C44  -  pv2)  =  0 


whose  three  solutions  give  the  magnitude  of  the  three  wave  speeds: 
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1/2 


(longi tudinal ) 
(transverse) 


(2.22) 


The  direction  of  propagation  of  these  three  waves  is  obtained  by 
evaluating  the  eigenvectors  a,  6,  y  from  Eq.  (2.21)  for  each  of  the 
speeds  v^ ,  ar|d  Vy  For  velocity  . v^,  we  find  a  =  1,  6  =  0,  and 

y  =  0.  A  comparison  of  a,  b,  y-with  I,  m,  n  shows  that  the  direction 
cosines  of  the  particle  displacement  are  identical  to  the  direction 
cosines  of  the  plane  wave  normal;  i.e.,  speed  v-j  corresponds  to  a 
pure  mode  longitudinal  wave.  Similar  calculations  for  a,  b,  y  for 
speeds  v^  and  v^  show  that  v 2  and  v^  correspond  to  pure  mode  transverse 
waves  and  that  | v^ 1  =  | v ^ 1 i  i.e.,  the  two  transverse  waves  propagate 
at  the  same  speed.  This  means  that  the  speed  of  a  transverse  wave 
propagating  along  the  [100]  direction  is  independent  of  polarization. 


For  the  [110]  plane,  the  direction  cosines  of  the  wave  normal 
are  a  =  1  / •/2",  m  =  ]/^2,  n  =  0.  For  this  direction, 

X11  =  A22  =  \  (C11  +  C44 ^ 

A12  A21  =  2  ^C12  +  C44 ^ 

(2.23) 

A13  =  A23  =  A31  A32  =  0 

A33  =  C44  ‘ 

Substitution  of  these  values  into  Eq.  (2.17)  and  simplification  gives: 

2^C1 1  "  C1 2 ^  "  pV^  f?(Cll  +  C12  +  2C44^  '  pv2]  (C44  “  pv2  =  0  .  (2.24) 
Solution  of  Eq.  (2.24)  gives  the  magnitude  of  the  three  wave  speeds: 

V1  =  {(C11  +  C12  +  2C44)/2p}1/2 

v2  =  {(C11  "  ci2)/2p}1/2  (2.25) 

rr  /  J/2 

v3  =  {C44/p} 

The  direction  of  propagation  of  these  three  waves  is  obtained 
by  evaluating  the  eigenvectors  a,  6,  and  y  from  (2.21)  for  each  of  the 
speeds  v^,  v^,  and  v^.  Note  that  in  this  case  v^  i  v^.  For  the  speed 
v -j  we  find  that  a  =  1  //2,  8  =  1//2,  y  =  0.  A  comparison  of  a,  8,  and 

y  with  £,  m,  and  n  shows  that  the  two  sets  are  identical,  so  v-|  is  a 

pure  mode  longitudinal  wave.  Similar  calculations  of  a,  6,  y  for  the 

speeds  and  v^,  and  comparison  with  £,  m,  and  n  shows  that  they  are 

pure  mode  transverse  waves. 


The  direction  cosines  of  the  wave  normal  in  the  [111]  direction 


are  2  =  m  =  n  =  1//3.  For  this  direction. 


Ml  "  x22  "  x33  "  2  ^11 


o  (C-n  +  ^^44^ 


(2.26) 


x12  =  x13  ”  X23  =  X21  =  X31  X32  =  3  ^C12  +  C44^  ' 


Using  Eqs.  (2.26)  and  (2.17),  the  magnitudes  of  the  three  wave  speeds 
are  given  by 


V1  = 


fcH  +  2C^  +  4C^] 


3p 


1/2 


v2  =  v3 


C11  ‘  C1 2  +  C44 


3p 


1/2 


(2.27) 


Again,  evaluating  the  eigenvectors  for  each  of  the  speeds  v-j ,  v^,  and 
v3  shows  that  for  v^  a  =  e  =  y  =  1//J. 

Comparison  of  a,  e,  and  y  with  1,  m,  and  n  shows  that  v^  is 
longitudinal.  Similar  comparisons  show  that  v2  and  v3  are  transverse 
and  equal  in  magnitude. 


D.  SUMMARY 

The  relationship  between  sound  velocity  and  C^'s  has  been 
derived  for  the  three  principal  directions  [100],  [110],  and  [111]  in 
a  cubic  crystal.  The  results  are  summarized  in  Table  2.1  where  it  is 
clear  that  in  the  [100]  and  [111]  directions  the  two  transverse  waves 
travel  at  the  same  speed,  but  in  the  [110]  direction  their  speeds  are 


Expressions  for  Velocities  and  K 


different.  For  the  purposes  of  this  thesis  one  is  interested  in  the 
relationship  between  C .  . 1 s  and  sound  velocity  for  longitudinal  waves 

*  J 

as  the  relationship  among  the  C^'s  for  longitudinal  waves  presented 
in  Table  2.1  are  identical  with  those  labelled  in  Table  1.1, 
p.  4.  To  emphasize  this  point  the  results  for  the  longitudinal 
waves  in  the  principal  direction  are  repeated  in  Table  2.2.  These 
are  the  relationships  to  be  analyzed. 


CHAPTER  III 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  configuration  used  in  the  measurement  of  the  velocity  of 
sound  and  the  related  Kg's  In  samples  of  germanium  is  shown  in 
Figure  3.1.  Electrical  pulses  modulated  at  30  MHz  are  generated  by  the 
gated  amplifier.  They  arrive  at  the  transducer  surface  through  the 
impedance  matching  network  and  cause  the  transducer  to  vibrate  and 
generate  a  30  MHz  pulse  of  ultrasonic  waves  which  propagate  through 
the  sample  and  are  received  by  the  capacitive  receiver  which  converts 
them  to  an  electrical  signal  which  is  amplified  and  displayed  on  an 
oscilloscope. 


A.  CAPACITIVE  RECEIVER  ASSEMBLY 

Details  of  the  capacitive  receiver  are  given  in  Figure  3.2.  The 
upper  end  of  the  receiver  assembly  consists  of  a  copper  electrode  1  cm 
in  diameter  surrounded  by  a  grounded  outer  assembly  that  is  Insulated 
from  the  electrode  by  a  Teflon  ring.  The  electrode  Is  spring-loaded 
to  make  contact  with  the  quartz  transducer  surface.  The  electrical 
signal  is  fed  to  the  quartz  transducer  through  the  spring. 

At  the  lower  end  of  the  assembly  is  an  electrode  of  1  cm 
diameter  placed  on  a  fused  silica  optical  flat.  A  grounded  concentric 
copper  ring  in  contact  with  the  optically  flat  end  of  the  sample 
provides  an  air  gap  between  the  electrode  and  the  sample,  a  separation 
of  the  order  of  10  microns.  When  the  sample  is  placed  in  position,  it 
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Detector 


measurements  using  the  pulse  overlap  technique 


touches  the  outer  ring.  The  optically  flat  sample  face  and  the  electrode 
surfaces  form  a  parallel  plate  capacitor  with  air  as  dielectric, 
as  shown  schematically  in  Figure  3.2.  The  capacitance  of  the  parallel 
plate  capacitor  may  be  evaluated  from 


where  A  is  the  area  of  the  plates  and  d  is  the  spacing  between  the 
plates  ( 5jj -1  Oy ) .  since  the  fringing  Is  negligible.  The  capacitance  is 
found  to  be  of  the  order  of  60-120  pF  for  a  spacing  of  d  =  5-10  microns 
and  for  detector  diameter  =  0.92  cm. 

B.  SAMPLES 

Measurement  of  the  velocity  of  sound  was  made  in  the  laboratory 
on  several  materials,  but  in  the  present  thesis  the  focus  is  on  samples 
of  germanium  single  crystals.  The  lattice  structure  of  germanium  is 
cubic  and  there  are  three  directions  along  which  pure  longitudinal 
ultrasonic  waves  propagate  In  cubic  lattices:  [100],  [110],  and  [111]. 
The  samples  allowed  measurements  to  be  made  along  these  pure  mode 
directions.  They  are  sketched  in  Figure  3.3  with  the  dimensions 
indicated.  The  samples  were  first  lapped  to  a  flatness  of  half  a  wave¬ 
length  of  light.  The  importance  of  making  the  samples  optically  flat 
results  from  the  plane  wave  assumption  used  in  the  theory  of  ultra¬ 
sonic  wave  propagation.  Since  we  used  high  frequency  (^30  MHz),  which 
corresponds  to  a  wavelength  of  0.17  mm  for  germanium,  even  a  small 
variation  in  the  surface  can  produce  a  significant  phase  shift  and 
affect  the  accuracy  of  the  experiment. 


After  the  samples  were  lapped,  they  were  washed  thoroughly  to 


remove  any  oil  that  might  have  Impeded  the  adherence  of  a  copper 
coating.  The  samples  were  placed  in  a  vacuum  of  10”^  Torr  and  allowed 

O 

to  outgas  for  one  hour.  Then  a  100  A  coating  of  copper  was  sputtered 
onto  the  surface  to  serve  as  an  electrode.  After  removal  of  the 
sample  from  the  vacuum  chamber,  a  quartz  transducer  was  attached  to 
one  of  the  optically  flat  conducting  surfaces  by  means  of  nonaqueous 
stopcock  grease,  then  the  sample  was  mounted  in  the  measurement 
assembly. 

The  output  of  the  cw  oscillator  (see  Figure  3.1)  was  applied  to 
the  gated  amplifier  whose  output  was  tested  with  an  oscilloscope  and 
applied  through  an  impedance  matching  network  across  the  quartz  trans¬ 
ducer  attached  to  the  sample  surface.  By  adjustment  of  the  impedance 
matching  network,  the  pulsed  oscillator-transducer  system  was  tuned 
for  maximum  power  transmitted  to  the  quartz  transducer  whose  resonant 
frequency  was  30  MHz. 

The  30  MHz  pulsed  ultrasonic  wave  generated  by  the  transducer 
travels  through  the  sample  and  is  reflected  between  the  sample  surfaces. 
When  it  reaches  the  lower  surface,  it  causes  the  sample  surface  to 
vibrate.  As  a  result,  the  capacitance  between  the  surface  and  electrode, 
as  shown  in  Figure  3.2,  changes  periodically.  This  produces  an  alternating 
current  which  is  amplified  and  observed  on  the  oscilloscope.  The 

gated  amplifier  is  capable  of  generating  a  second  pulse  train  (delayed 
in  time  with  respect  to  the  first  pulse  train).  The  two  pulsed  ultra¬ 
sonic  waves  interfere  and  give  rise  to  interference  maxima  and  minima. 
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For  a  minimum,  the  path  difference  between  the  two  pulses  must  contain 
an  odd  number  of  half-wavelengths;  viz, 

Path  difference  =  (m  +  j)  x 

where  m  =  0,  1,  2,  _  In  taking  data  an  initial  frequency  in  the 

range  of  the  resonance  frequency  was  chosen.  Then  a  number  of  minima 
was  counted  as  the  frequency  was  increased,  and  the  final  frequency 
noted. 

If  the  nth  echo  of  one  pulse  train  overlapped  the  mth  echo  of 
the  second  pulse  train,  the  difference  m  -  n  =  s  is  the  quantity  used 
in  calculating  the  velocity. 

Finally,  the  length  L  of  the  sample  was  measured  with  an  accuracy 
-4 

of  10  cm.  The  velocity  of  the  ultrasound  in  the  sample  material  was 
then  determined  using  the  expression: 

v  =  ^S...AF 
AQ 

where 

L  =  the  length  of  the  sample; 

s  =  m  -  n  (n  is  the  echo  number  of  some  Initial  pulse  and  m  is 
the  echo  number  of  some  other  pulse  delayed  in  time  with 
respect  to  the  first  pulse); 
aF  =  frequency  difference; 
aQ  =  the  number  of  interference  minima. 


CHAPTER  IV 


EXPERIMENTAL  RESULTS  AND  DATA  ANALYSIS 

The  error  propagation  in  determining  the  K^'s  from  expressions 

for  C..  can  be  evaluated  directly  from  the  data;  however,  in  order  to 

evaluate  the  relative  magnitude  of  the  propagated  error  and  the  error 

resulting  from  direct  measurement  of  the  l^'s  it  was  necessary  to  set 

up  apparatus  and  measure  the  velocity  of  longitudinal  waves  in 

crystalline  samples.  The  calculations  given  in  Chapter  II  show  that 

the  expressions  for  the  l^'s  are  related  to  the  measured  velocities 
? 

by  =  pv  .  This  chapter  describes  the  means  by  which  experimental 
values  of  the  longitudinal  wave  velocities  were  obtained. 

A.  EXPERIMENTAL  DATA 

The  length  and  the  density  of  the  samples  used  in  the  present 
experiment,  viz.  Ge(100),  Ge(llO),  and  Ge(lll),  are  given  in  Table  4.1. 
Three  sets  of  data,  using  30  MHz  transducers  were  taken  for  the  measure¬ 
ment  of  velocity  of  sound  in  the  samples  of  germanium.  For  each  set  of 
data,  a  measurement  of  frequency  F  was  taken  for  several  values  of  the 
maxima  Q.  The  data  are  presented  in  Tables  4.2,  4.3,  and  4.4.  Plots 
of  F  versus  Q  are  shown  for  each  of  the  samples  in  Figures  4.1,  4.2, 
and  4.3.  These  figures,  usually  called  scatter  diagrams,  allow  one  to 
observe  directly  the  scatter  of  the  data  resulting  from  statistical 
variations  in  the  magnitudes  of  the  measured  quantities.  Under  ideal 
circumstances  these  variations  would  vanish  and  the  data  in  each  case 


for  Ge[100] 
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Figure  4.2.  Scatter  diagram  of  F  vs.  Q  for  Ge[110]. 
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Figure  4.3.  Scatter  diagram  of  F  vs.  Q  for  Ge[lll]. 


would  define  a  single  straight  line  whose  slope  in  combination  with  L 
and  s  would  give  the  exact  value  of  the  velocity  in  the  corresponding 
direction.  The  purpose  of  the  next  section  is  to  define  the  way  one 
handles  the  actual  data  to  arrive  at  the  most  probable  straight  line, 
one  that  is  a  very  good  approximation  to  the  exact  one. 

B.  DATA  ANALYSIS 

To  calculate  the  velocity  of  sound  in  the  germanium  crystal 
for  various  crystalline  modes  one  uses  the  expression 

»=2  Lsf£  (4.,) 

and  the  observations  on  L,  s  (m  and  n),  F,  and  Q  presented  above. 

The  observations  on  F  and  Q  presented  in  the  scatter  diagrams 
(Figures  4.1,  4.2,  4.3)  suggest  that  for  any  interval  aQ  the  corres¬ 
ponding  interval  aF  is  not  uniform  for  all  values  of  Q.  For  example, 

if  AQ  =  Q2  -  Q.|  and  aF  =  F2  -  F-| ,  then  for  a  different  aQ  = 

(=  Q2  -  Q-jh  F^  -  F^  can  be  different  from  AF  =  -  F-| .  This  is 

because  the  data  on  F  and  Q  in  practice  do  not  have  an  exact  functional 

relationship  as  a  result  of  systematic  and  random  errors  in  the  experi¬ 

ment.  It  therefore  is  important  statistically  to  estimate  the  best 
value  of  aF/aQ  from  the  data  in  order  to  predict  the  most  probable 
value  of  the  velocity  v.  This  is  the  purpose  of  the  present  section 
which  is  divided  into  five  parts. 

The  first  part  deals  with  the  assumptions,  relevant  formulae, 
and  methodology  employed  in  the  estimation  of  aF/aQ.  In  the  second 
part,  ordinary  least  squares  estimates  of  the  slope  6  =  aF/aQ  are 


presented  in  the  standard  format  and  also  in  a  table.  Computations 
of  velocity  are  then  made  and  the  results  are  shown  in  tabular  form. 

In  the  third  part,  the  experimental  results  on  velocity  in  the 


germanium  crystals  are  compared  with  those  obtained  by  McSkimin  (1963). 
The  propagation  of  errors  in  the  present  experiment  and  propagation  of 
error  in  calculating  velocity  using  standard  values  of  second-order 
elastic  constants  are  discussed  in  the  fourth  part.  Lastly,  in  the 
fifth  part  are  presented  the  conclusions. 

1 .  Assumptions  and  Methodology 

The  method  of  ordinary  least  squares  is  used  to  estimate  the 
most  probable  value  of  the  slope  B  =  AF/AQ  from  the  data.  But  in 
order  to  do  so  it  is  essential  to  make  the  (unverifiable)  assumption 
that  the  statistical  distribution  of  F  does  not  change  from  one  set 
of  observations  to  another.  The  evaluation  of  the  slope  B  requires 
a  specification  of  the  relationship  F  =  f(Q)  and  the  use  of  statisti¬ 
cal  regression  analysis  [Kleinbaum  and  Kupper  (1978);  Johnston  (1984)]. 

In  what  follows  a  possible  relationship  between  F  and  Q  is  assumed  and 
a  methodology  for  statistical  estimation  of  B  is  developed. 

a.  Functional  relationship  between  F  and  Q.  The  scatter  diagrams 
in  Figures  4.1,  4.2,  and  4.3  show  that  in  observations  on  each  sample 
the  value  of  F  at  a  certain  value  of  Q  is  not  consistent.  This  is  due 
to  the  fact  that  random  errors  impart  bias  to  the  observation  on  F. 

In  fact,  the  distribution  of  F  at  the  same  value  of  Q  depends  upon  the 
statistical  distribution  of  the  random  error.  Thus,  there  is  a  whole 
probability  distribution  of  values  of  F  for  each  value  of  Q.  The 


scatter  diagrams  further  suggest  that  the  variation  in  the  mean  of 
F  values  at  some  Q  is  approximately  linear  with  Q.  Thus,  a  stochastic 
linear  relationship  between  F  and  Q  can  be  assumed.  This  stochastic 
linear  relationship  would  become  deterministic  if  the  variance  of  F 
were  zero;  i.e.,  if  there  were  no  random  errors  in  observations  of  F. 
The  stochastic  relationship  between  F  and  Q  for  n  observations  is: 

F.j  =  a  +  6Qi  +  i  =  1,  2,  ...  n  (4.2) 

where  the  observed  frequency  F^  is  taken  to  be  the  dependent  variable, 
the  number  of  maxima  is  the  independent  variable,  and  t  is  the 
stochastic  disturbance,  a  and  g  are  the  regression  parameters.  The 
values  of  F..  and  are  observable  but  those  of  are  not.  The 
purpose  of  the  term  is  to  characterize  the  discrepancies  that  emerge 
between  the  actual  observed  values  of  F  and  the  values  that  would  have 
been  given  by  an  exact  functional  relationship  of  the  form: 

F  =  a  +  6Q  .  (4.3) 

The  fact  that  the  magnitude  of  e..  cannot  be  determined  exactly  means 
that  the  value  of  F  can  never  be  forecast  with  certainty,  i.e.,  with 
probability  1.  The  uncertainty  concerning  F..  arises  due  to  the 
presence  of  the  stochastic  disturbances  which,  being  random,  imparts 
randomness  to  F. .  The  randomness  in  the  term  may  be  on  account  of 
a  variety  of  factors  which  may  or  may  not  all  be  quantifiable.  Among 
those  factors  are  random  apparatus  or  human  ones  as  well  as  systematic 
errors  in  measurement  of  F,  L,  and  Q.  The  net  effect  of  all  such 


factors,  then,  is  summarized  by  a  single  stochastic  variable  .  The 
probability  distribution  of  F  and  its  properties  are  then  determined 
by  the  values  of  Q  and  by  the  probability  distribution  of  e.  Thus,  the 
full  description  of  the  model  in  Equation  (4.2)  also  calls  for  the  full 
specification  of  the  probability  distribution  of  the  error  .  We 
make  the  following  statistical  assumptions: 

1.  is  normally  distributed  V-  i . 

2.  e.j  has  zero  mean,  i  .e. ,  <e1->  s  Olf  i. 

2  2 

3.  e.j  is  homoskedastlc,  i.e.,  <e^  >  =  a  .  This  means  that 

2 

every  disturbance  has  the  same  variance  o  for  all  observa¬ 
tions  whose  value  is  unknown.  This  assumption  rules  out,  for 
example,  the  possibility  that  the  disturbance  could  be  greater 
for  higher  values  of  Q  than  for  lower  values  of  Q. 

4.  e.  is  nonautoregressive;  i.e.,  <e.e.>  =  0  if  i  f  j.  This 

1  1  J  f  0  i  =  j 

assumption  implies  that  the  expected  value  of  F  at  any  time 

in  an  experiment  will  be  different  from  the  expected  value 

of  F  in  the  same  experiment  at  a  different  time. 

5.  Q.j  is  a  nonstochastic  variable  with  values  fixed  in  repeated 

sets  of  observations  such  that  for  any  sample  size 

1  n  _ 2  _ 

—  I  (Q.  -  Q)  is  a  finite  number  where  Q  is  the  mean  value 
n  1=1  1 
of  Q. 

With  these  assumptions,  we  can  find  the  properties  of  the 
probability  distribution  of  F^  for  all  i.  The  mean  of  F^  is  obtained  by 
taking  the  mathematical  expectation  value  of  both  sides  of  Eq.  (4.2): 


f a . -a  -a  a •  .  \  *\  a  -a  a a  .  a «.  -a  a- a  •/ 


Since  the  expected  value  of  ^  vanishes  ;i .e. ,  <e^>  =  0,  then 


<F^>  =  <3!  +  eQ^>  . 


(4.4b) 


The  variance  of  F.  is 


Var(F.)  =  <[Fi  -  <F.>]S 

=  "[(a  +  6Qi  +  ) 

2  2 

=  <e.j  >  =  a 


-  (» -  sQi)r> 


(4.5) 


Equation  (4.4b),  which  gives  the  mean  value  of  F  for  each  value 
of  Q,  is  known  as  the  population  regression  line.  The  slope  g  measures 
the  change  in  the  mean  value  of  F  corresponding  to  a  unit  change  in  the 
value  of  Q. 

Estimation  of  the  values  of  a  and  8  gives  the  sample  regression 
line  that  serves  as  an  estimate  of  the  population  regression  line.  If 
a  and  8  are  estimated  by  a  and  S,  respectively,  then  the  sample 
regression  line  is 


^  -  a  +  6Q1  (4.6) 

where  is  the  estimate  of  F.  or  the  fitted  value  of  F^ .  The  observed 
values  of  F^  does  not  necessarily  lie  exactly  on  the  sample  regression 
line  so  that  the  value  of  F.  and  F^  in  general  are  different.  This 


difference  is  called  the  residual  and  is  denoted  by  e^ .  Thus  we  dis¬ 
tinguish  between  the  following: 


( 


F..  =  a  +  8Q..  +  (population)  , 

Fi  =  a  +  8Q^  +  ei  (sample)  . 


(4.7) 


In  general,  the  residual  e^  is  different  from  because  a  and  6  are 
different  from  the  true  values  of  a  and  8.  In  fact,  e.  is  the  estimate 
of  the  disturbance  e^.  Figure  4.4  is  a  schematic  representation  of 
this  distinction.  The  slope  6  is  obtained  from  Eq.  (4.6): 


^1.  a£.  e 

dQ  aQ  b 


(4.8) 


b.  Evaluation  of  intercept  a  and  slope  B.  a  and  8  are 
estimated  using  ordinary  least  squares.  If  each  residual  e^  is 
squared,  negative  signs  disappear,  and  the  sum  of  squared  residuals 
is  a  nonnegative  quantity.  In  using  the  least  square  principle  one 
selects  a  and  8  for  minimum  ie/.  First,  one  evaluates 


>.  =  F-  -  F.  =  F.  -  (a  +  80.)  , 

l  l  i  l  v  ’ 


ze.2  =  E(Fi  -  a  -  §Qi )2  . 


The  necessary  conditions  for  a  stationary  minimum  are: 


=  0  , 


=  0  . 


(4.9) 


v  r  <V-  *  -  V  ‘  / 


SVN  Ev. 


It  follows  that: 


IF.  =  na  +  6  Z  Qi 


where  n  is  the  number  of  data  points.  Since  nQ^  =  Z  ,  the  product 


II  II  A  II  q 

l  F-j  Q-j  ■  a  I  Qi  +  B  I  Qi  • 

i =1  1  1  i=l  1  1=1  1 


In  matrix  form  this  becomes 


'a 

s 

t-' 

l 

n  zQ,  "  f zF. 


'i  lEFi«iJ 


Evaluating  the  elements  of  the  matrix,  one  obtains 


“  -  H- £  Fi  -  5  l  I  . 


(4.10a) 


s(Qi  -  Q)(F1  -  F) 


i(Qi  -  Q)2 


(4.10b) 


The  quantities  a  and  8  are,  respectively,  the  intercept  and  the  slope 
of  the  regression  line.  The  slope  8  is  used  to  calculate  the  mean 
value  of  the  velocity  reported  in  this  thesis. 


2.  Use  of  Data  in  Ordinary  Least  Squares  Estimates 


of  e  and  Velocities 


a.  Computer  program.  A  computer  program  for  evaluating  the 
regression  line  is  available  at  The  University  of  Tennessee  Computer 
Center.  This  program,  SYSREG,  presented  in  Appendix  A,  was  used  to 


stm.  tC M.  mjf* AjV--,  V-/ .O  .t  ft  A. Jk*, 


evaluate  the  sample  regression  lines  as  In  Equation  (4.7).  Values  of 
a  and  6  evaluated  from  this  regression  line  customarily  are  presented 
as  follows: 


F.  =  a  +  |Q1  R2  =  ...  (4.11) 

(SO  (SO 

a  D 

2  2 
where  Sa  is  the  estimate  of  variance  of  a  and  Sg  is  the  estimate  of 

„  2 

variance  of  8.  R  is  the  coefficient  of  determination  which  is  a 

measure  of  "goodness  of  fit";  i.e.,  how  well  the  sample  regression  line 

2 

fits  the  observations.  R  indicates  the  proportion  of  variation  of  F 
that  can  be  attributed  to  the  variation  of  Q.  It  is  evaluated  from 

R2  _  Regression  Sum  of  Squares 
Total  Sum  of  Squares 

.  1  * 
l  (q  - F) 

2  2  2  • 

R  takes  on  the  values:  0  <_  R  <_  1 .  A  zero  value  of  R  indicates 

the  poorest  and  a  unit  value  the  best  fit  that  can  be  attained. 

b.  Velocities  in  germanium.  The  estimated  regression  lines 
obtained  by  applying  the  ordinary  least  squares  estimate  to  data  on 
the  samples  Ge(lll),  Ge(110),  and  Ge(100)  are  shown  in  standard  form 
in  Table  4.5.  Note  that  the  value  of  R2  in  all  the  three  orientations 
is  close  to  unity  which  indicates  that  the  data  are  well  fit  by  the 
regression  lines.  The  results  are  repeated  in  more  detail  in  Table  4.6, 


Table  4.5.  Estimated  Regression  Lines  for  Ge[100],  Ge[110],  and 
Ge[lll] 


Sample 

Ge(100) 


Ge (110) 

F,  =  24.916158  +  0.1118  Q.  +  e.,  R2  =  0.9996 

1  11 

(0.0134)  (0.0003) 

Ge  ( 11 1 ) 

F.  =  25.7741  +  0.05694  Q.  +  e-  ,  R2  =  0.9992 
(0.0134)  (0.000156) 

_ Estimated  Regression  Line _ 

Fi  =  27.8650  +  0.0660  Q1  +  e.  ,  R2  =  0.9986 

(0.0133  (0.0002) 


where  a  is  the  intercept  of  the  estimated  regression  line,  8  is  the 

slope,  S-  and  S~  are  the  standard  deviations  in  the  measurements  of 

a  8 

a  and  8. 


3.  Comparison  of  Results  with  Reference  Values 


The  comparison  of  the  results  of  the  present  measurements 
presented  in  Table  4.6  with  those  of  McSkimin  (also  presented  in 
Table  4.6)  is  especially  informative.  The  present  data  present  the 
scatter  resulting  from  all  sources  of  error,  both  definable  and 
undefinable.  Among  the  definable  sources  of  error  are  systematic 
errors  resulting  from  measurement  of  sample  length,  resonant  fre¬ 
quency,  and  density.  Random  errors  from  repetition  of  these  measure¬ 
ments  also  enter. 


a.  Evaluation  of  scatter  of  data  around  the  mean.  The  results 
on  velocity  reported  in  Table  4.6  give  the  deviation  of  velocity  values 
from  the  mean  for  a  95%  probability  level.  The  specification  of  the 
probability  level  gives  a  more  complete  picture  of  the  effect  of 
random  errors  on  the  data  than  usually  is  given.  In  Table  4.6  our 
estimated  velocity  in  Ge[  1 11]  falls  in  the  interval  ±  0.53%  of  the  mean 
value  of  velocity  with  a  95%  probability.  McSkimin  reported  a  mean 
square  deviation  of  ±  0.02%.  It  is  apparent  that  our  mean  square 
deviation  is  larger  than  McSkimin's,  but  the  significance  of  this 
statement  is  somewhat  difficult  to  evaluate  without  information  about 
his  probability  level.  Further,  whether  the  difference  results  from 
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our  consideration  of  all  sources  of  random  error  while  McSkimin 


considered  only  errors  resulting  from  diffraction  is  equally 


uncertain.  Comparison  of  the  remaining  data  in  Table  4.6  results  in 


similar  conclusions. 


b.  Comparison  of  our  mean  value  of  velocity  with  values 


liven  by  McSkimin.  It  remains  to  compare  the  mean  value  of  velocity 


in  this  experiment  with  the  values  of  velocity  given  by  McSkimin. 


Agreement  between  the  two  sets  of  values  adds  credibility  to  the 


results.  We  let  the  three  velocities  in  a  germanium  crystal  as 


measured  by  McSkimin  be  v^-],  v^g  and  v-jqq-  Somewhat  at  random, 
we  chose  the  velocity  in  the  [111]  direction  for  our  discussion  of 
the  comparison.  (The  comparison  for  the  [110]  and  [100]  directions 


is  obtained  by  using  velocity  values  appropriate  to  those  directions.) 


Corresponding  to  McSkimin's  value  we  evaluate  the  slope 


(4.12) 


To  test  whether  the  velocity  estimated  in  this  study  is  equal  to  that 


of  McSkimin,  we  have  to  test  whether  our  B  is  equal  to  b-jii-  There¬ 


fore,  we  set  the  null  hypotheses,* 


Ho:  8  s  6nr 


*A  null  hypothesis  is  a  proposition  which  is  considered  valid 
unless  evidence  throws  serious  doubt  on  it. 


-ts 


Equivalently, 


H0:  v  =  V111- 


The  T-statistic  for  the  [111]  direction  is  the  function 


The  distribution  of  the  T-statistic  about  zero,  known  as  the  t-distri- 
bution,  has  (n-2)  degrees  of  freedom.  In  estimating  F  one  has  two 
unknowns  a  and  6  whose  presence  reduces  the  numDer  of  degrees  of  free¬ 
dom  from  n,  the  total  number  of  observations,  to  n-2.  Let  the  signifi¬ 
cance  level  (unity  minus  the  probability)  be  designated  A.  Then,  for 
a  specified  significance  level  A,  the  T-statistic  lies  between  the  lower 
limit  and  the  upper  limit  of  a  critical  region  as  follows: 


6  -  3 


-  t 


111 


<  t. 


n-2, A/2  -  -  -  n-2, A/2 

B 


(4.13) 


The  critical  region  is  shown  in  Figure  4.5  with  the  limits  specified. 
Without  the  limits  the  t-distribution  would  vary  from  -°°  to  <*>.  In  that 


case  the  cumulative  probability  would  be  unity;  i.e.,  the  area  under 

fe  -  s, . 

is  unity.  For  a  significance  level  A 


the  density  function  f 


’111 


B 


the  area  under  the  distribution  function  is  the  probability 


<■  n-2 


.a/2 


P  = 


-t 


n-2, A/2 


B  -  6 


111 


B  -  8 


111 


S6 


s 


which  usually  is  designated  as  follows: 
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8  -  0, 


P  =  -tn-2,x/2  -  S'  -  tn-2,X/2  =  1  ' 

P 


(4.14) 


The  subscripts  on  the  limits  include  x/2  because  an  area  of  x/2  is 


found  on  each  tail  of  the  t-distribution  between  t  0  and  «. 

n-2 ,  X/ 2 


These  areas  can  be  read  from  t-distribution  tables  In  statistics  text¬ 


books  [Kleinbaum  and  Kupper  (1978)].  For  data  in  this  thesis  a  value 


X  =  1%  is  chosen,  meaning  that  the  area  of  the  two  tails  is  0.5%  each. 


The  corresponding  critical  region  and  the  T-statistic  for  the  three 


crystalline  modes  are  evaluated  and  are  reported  in  Table  4.7. 


From  the  inequality  (4.13)  we  have 


Blll  "  S8  tn-2,X/2  -  B  -  6111  +  SB  tn-2,X/2 


(4.15) 


The  above  inequality  gives  the  lower  and  upper  limits  of  fj  for  the 


chosen  significance  level  X.  The  corresponding  lower  and  upper 


limits  of  the  estimated  velocity  v  are  found  from  (4.15)  by  multiply¬ 


ing  the  inequality  by  a  factor  2Ls: 


vlll  "  2Ls  *  tn-2,x/2  *  S§  1  V  1  vlll  +  2Ls  ’  tn-2 ,x/2  *  Se  ,(4J6) 


terms  of  a  percentage  of  we  have: 


L  tn-;,>/2  s6x100v  <v<» 

“ill  - 111 


( t  ?  •  S') 

+  2Ls  n~2».x/i - 8- 


x  100  v. 


(4.17a) 


■aw-# -ys&te >:■ 


SC 
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vm  -  Avni  i  v  -  vm  +  Avm 


(4.17b) 


where 


=  [2LS  t-n-2-^^-2  x  II 


Similar  results  are  evaluated  for  the  other  two  crystalline  modes. 

The  results  are  presented  in  Table  4.7,  column  7.  The  associated 
probabilities  are  calculated  from  (4.14)  and  the  results  are  presented 
in  Table  4.7,  column  8. 

From  the  results  reported  in  Table  4.7  we  find  that  for  x  =  1% 
or  for  a  probability  level  99%  the  T-statistic  does  not  exceed  the 
critical  region  for  all  the  three  crystalline  modes.  As  long  as  the 
T-statistic  is  within  the  critical  region  .for  a  specified  significance 
level,  we  do  not  reject  the  null  hypothesis  and  conclude  that  there  is 
not  enough  evidence  to  suggest  that  b  is  different  from  6^,  or  that 

the  equivalently  estimated  velocity  v  is  different  from  McSkimin's 

velocity  In  other  words,  the  deviation  of  v  around  v^^  will 

be  confined  to  upper  and  lower  limits  v^  ±  Av^  for  99%  of  the 

time.  Only  1%  of  the  time  the  deviation  of  v  from  v^^  will  be  large 

enough  so  as  to  fall  outside  the  limits  v^^  ±  Av^^.  Similar  or 
better  results  are  obtained  for  the  other  two  crystalline  modes. 

Thus  we  may  conclude  that  the  estimates  of  velocity  in  our  experi¬ 
ment  are  in  agreement  of  those  of  McSkimin.  The  results  add  credibility 
to  our  experiment  and  results  as  well  as  McSkimin's  results. 


■Wv 
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4.  Error  Propagation 

In  this  section  we  compare  the  relative  error  in  estimation 

of  velocity  in  the  present  study  with  the  error  obtained  by 

recombining  values  of  elastic  constants  measured  by  others  and  using 

2 

the  expression  =  pv  . 

a.  Error  propagation  in  the  present  experiment.  To  derive  a 
formula  for  the  percentage  error  in  the  velocity  one  uses  the 
expression  for  the  regression  line  (Eq.  4.7). 

Fi  =  a  +  6Qi  +  e1  ,  (4.18) 

where  i  =  1,  2,  ...  n.  After  applying  a  least  squares  analysis  the 
fitted  value  of  is 

Fi  =  a  +  6Qi  .  (4J9) 

The  mean  values  of  F^  and  also  fit  the  regression  line.  Therefore, 
we  have 

r  =  s  ♦  ar  .  (4-20> 

where  T  and  Q  are  the  mean  values  of  the  observations  taken  in  the 
experiment  for  a  particular  crystalline  mode.  These  values  are  given 


as  follows: 


From  Eq.  (4.20), 


B  =  - 


F  -  a 


Taking  the  logarithm  of  both  sides  and  evaluating  the  partial 


derivatives: 


i!  =  3(F  -  a)  _  3C[ 
B  (F  -  a)  Q 


(4.21) 


Evaluation  of  the  relative  errors  requires  consideration  of  the  absolute 


value  of  each  term: 


66  _  6 ( F  -  a )  +  6Q 
8  (F  -  i)  Q 


(4.22) 


Since  v  =  2Lsb,  the  relative  error  in  velocity  is:  —  =  —  as  L  and 


s  are  constants.  Therefore, 


6F  -  S~ 

6v_ _ i  +  iH. 

V  "  F  -  a  Q 


(4.23) 


F  -  a 


kSjl 


V  >  V 


i 


I 


I 


■ 

>, 

V 


i 
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where  the  identity  S-  =  S  has  been  used.  This  formula  is  used  to 

a  a 

evaluate  errors  in  the  measurement  of  velocity  due  to  random  causes  for 
all  three  samples  and  the  results  are  presented  in  Table  4.8. 

The  error  largely  depends  upon  the  value  of  the  mean  F.  The  larger  the 
value  of  F,  the  smaller  the  percentage  error.  The  error  in  Ge  (111)  is 
much  smaller  than  the  error  in  Ge  (100)  because  it  was  possible  to 
obtain  a  larger  number  of  measures  in  the  Ge  (111)  sample,  and  hence 
to  have  a  larger  value  of  F.  The  experimenter  does  not  always  have  total 
control  over  the  value  of  F,  but  should  always  seek  to  obtain  the  largest 
value  possible. 


b.  Error  propagation  in  use  of  reference  values  of  C^.  The 
percentage  error  found  in  the  preceding  section  now  is  compared  with 
the  error  propagated  when  one  takes  the  values  of  the  second-order 
constants  from  reference  sources  (McSkimin,  1963)  and  calculates  the 
velocity  of  longitudinal  waves  in  the  three  principal  directions  using 
the  following  formulae: 


ioo  =  cn/p 

(4.24a) 

110  =  (Cll  +  C12  +  2C44^2p 

(4.24b) 

111  ”  ^11  +  2^12  +  ^^44^^p 

(4.24c) 

in  which  the  combinations  of  second-order  elastic  constants  are 
recognized  as  being  the  same  as  the  l^'s  listed  in  Table  1.1,  p.  4). 
The  relative  errors  in  the  velocities  are: 


I 
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Table  4.8.  Propagated  Error  (6v/v)  in  Present  Experiment 


Sample 

Propagated  Error  (5v/v) 

Ge [ 1 00] 

0.46% 

Ge[110] 

0.33% 

Ge[ 1 11] 

0.29% 

I 

& 

[5 

K 

E 

i 

V. 

I 

E 

I 


i 

& 

V- 


6v100 
V1 00 


5C 


11 


'11 


1_  5e^ 
2  p 


6vno 

vno 


i_ 

2 


fiC-ji  +  <5  C 1 2  +  26C^ 


C11  +  C12  +  2C44 


X  Xe. 

2.  p 
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(4.25a) 

(4.25b) 


111 

11 


6C-|i  +  2  6C 1 2  +  4SC^ 
C11  +  2C12  +  4C44 


+  1  <5p 

2  p  ' 


(4.25c) 


Using  these  equations  and  the  values  for  the  errors  given  in  the 
reference  and  assuming  an  error  in  the  measurement  of  p  as  0.1%,  the 
propagated  errors  are  calculated  and  are  given  in  Table  4.9. 

5.  Correction  of  Systematic  Error 

Close  examination  of  the  data  in  Tables  4.2,  4.3,  and  4.4 
(pp.  28,  29,  30)  presented  in  the  scatter  diagrams  in  Figures  4.1, 

4.2,  and  4.3  (pp.  31,  32,  33)  reveals  systematic  errors  in  some  of 
the  data.  Some  of  the  numbers  appear  to  deviate  by  more  than  one 
standard  deviation  resulting  from  either  an  extra  count  or  a  missed 
count  in  the  data  taking.  Such  errors  can  be  corrected  by  standard 
techniques.  The  data  were  subjected  to  analysis  by  a  computer  program 
to  make  such  corrections.  The  results  are  given  in  Tables  4.10  and  4.11 
which  can  be  compared  directly  with  Tables  4.6  and  4.7  (pp.  44,  50), 
respectively.  Comparison  reveals  that  both  the  slopes  6  and  the  standard 
deviations  Ss  were  improved.  The  slopes  were  used  to  calculate  present 

D 

experiment  values  of  the  velocities  v  which  are  observed  to  be  in  better 
agreement  with  the  measured  data  of  McSkimin  listed  in  Table  4.10,  and 
the  deviations  from  the  mean  velocity  were  reduced.  In  Table  4.11  it 


l 


Cu  =  12.8528  x  10"  ±  0.04% 
C22  -  4.8259  x  10"  ±  0.04% 
C44  =  6.67966  x  10"  ±  0.04% 


Comparison  of  Corrected  Velocities  in  Present  Experiment  with  Those  of  McSkimin 


is  found  that  the  T  statistic  still  lies  within  the  critical  region 
for  the  [100]  and  [110]  values;  however,  the  [111]  T  statistic  now  lies 
outside  the  critical  region.  This  probably  results  from  the  fact  that, 
as  indicated,  the  [111]  value  was  obtained  from  McSkimin's  data  by 
adding  certain  numbers  to  give  a  propagated  error  which  has  not  been 
accounted  for  in  the  analysis.  An  additional  possibility  is  that  the 
present  value  of  Ss  of  only  0.000053  results  in  an  anomalously  large 

P 

T  statistic.  Finally,  the  variation  Av  of  the  present  data  from  those 
of  McSkimin  in  all  three  cases  is  decreased  by  the  correction  of  the 
systematic  error. 

The  results  presented  in  Tables  4.10  and  4.11  also  have  been 
calculated  by  including  all  of  the  significant  figures  in  the  lengths 
given  in  Table  4.1  (p.  27).  Tables  4.6  and  4.7  (pp.  44,  50)  were  calcu¬ 
lated  by  rounding  off  the  lengths  to  three  significant  figures.  Although 
the  change  in  the  velocities  resulting  from  the  more  accurate  value  of 
length  is  not  great,  it  was  detectable  in  the  fourth  significant  figure 
in  the  velocity.  Hence,  the  correction  is  justified. 


CHAPTER  V 


SUMMARY 

The  analysis  in  this  thesis  shows  that  the  question  originally 
posed  does  not  have  a  unique  answer  for  all  samples  under  all  con¬ 
ditions.  One  cannot  decide  a  priori  whether  reference  values  of  C-. 

^  J 

should  be  used  or  whether  one  should  measure  the  C..  each  time 

1 0 

he  measures  the  C^.  The  analysis  given,  however,  tends  to  support 
the  position  that  on  those  occasions  one  has  data  as  accurate  as  those 
of  McSkimin  his  accuracy  is  greatest  if  he  uses  them  rather  than 
remeasuring  each  sample.  If  such  accurate  data  are  not  available,  one 
has  no  choice. 
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APPENDIX  A 


THEORY  OF  PULSE  SUPERPOSITION  TECHNIQUE 


The  equation  of  motion  for  a  progressive  wave  propagating  in  a 


medium  is  given  by 


u  =  Aei{kx-wt) 


(A-l) 


where  A  is  the  amplitude  of  the  wave. 

Let  us  consider  two  such  progressive  waves  given  by 


u,  =  A,e 


i  (kx^  -wt-j ) 


u2  =  A2e 


i  (kx2-wt2) 


(A-2) 


Consider  the  sample  of  length  i  as  shown  in  Figure  A.l.  The  electri¬ 
cal  signal  applied  to  the  transducer  bonded  to  the  sample  by  means  of 
stopcock  grease  causes  the  transducer  to  emit  an  ultrasonic  wave  which 
travels  through  the  sample  and  is  detected  at  the  opposite  end  of  the 
sample  by  the  capacitive  receiver.  The  signal  is  displayed  on  the 
oscilloscope  screen.  In  this  process  the  first  pulse  is 
detected  when  the  ultrasonic  wave  reaches  the  end  of  the  sample  at 
the  capacitive  receiver.  After  this,  the  wave  undergoes  reflection, 
returns  and  once  again  is  reflected.  The  second  pulse  thus  seen  on 
the  oscilloscope  screen  accounts  for  the  signal  that  is  detected 
after  the  ultrasonic  wave  has  undergone  two  reflections. 


Since  the  expression  for  the  velocity  of  sound  is  derived  upon 
the  basic  assumption  that  the  wave  has  undergone  two  reflections 
before  it  is  detected,  the  velocity  formula  derived  below  is  valid  for 
all  pulses  except  the  first  pulse  which  is  detected  before  two 
reflections.  Therefore,  for  the  first  echo  of  some  initial  pulse  which 
corresponds  to  the  second  pulse  that  is  seen  on  the  screen  as  shown  in 
Figure  A-2,  one  can  write 

x]  =  (x  +  2i)  +  x  ’  (A-3) 

The  pulse  travels  the  distance  £,  is  reflected  at  the  sample  surface 
and  again  travels  a  distance  £,  whereupon  it  Is  reflected  again  at  the 
boundary.  After  the  second  reflection,  if  it  travels  a  distance  x,  as 
shown  in  Figure  A-3,  then  the  total  distance  traveled  is  (2'  +  x).  In 
addition,  since  the  wave  undergoes  two  reflections,  it  also  undergoes 
two  phase  changes  which  need  to  be  accounted  for.  The  distance  corres¬ 
ponding  to  one  phase  change  is  £,  where  $  is  the  phase  change  upon 
reflection.  Since  two  phase  changes  are  involved,  the  corresponding 
distance  Is  given  by  x*  Therefore,  the  total  distance  traveled  by 
the  wave  before  the  first  echo  is  seen  is  given  by  the  sum 

X-|  =  x  +  h  +  . 

If  there  are  n  such  echoes,  the  expression  takes  the  general  form 


for  the  nth  echo.  Similarly,  for  some  other  Initial  pulse  (delayed 
in  time  with  respect  to  the  first  pulse)  for  the  mth  echo,  we  have 


x2  =  x'  +  2m  +  . 


For  overlap  of  echoes  we  require 


i(kx1-wt1)  i(kxp-a)t?) 
u  =  u1  +  u2  =  A^e  +  Ape  /" 


(A- 


Since 


we  have 


u  =  Agi  (kx-aitjj-gi  (2nk£+2n<j>)  +  0i  (2mk£+2m4> ) -| 


Since  the  delayed  pulse  started  later  than  the  first  one,  we  have 
n  >  m,  so  that  one  can  write 


Substituting  for  n  in  terms  of  m  and  s,  we  have 


u  =  Ae1(kx-wt)|-ei{2(m+s)kt +2(m+s)4>} +gi(2mkt +2m0)-| 


u  =  Aei(kx-u)t)ei(2mki  +  2m4))j-ei(2sk£+2s<{.)  +  ^  ^A_?^ 


For  destructive  Interference,  we  have  the  condition 


2sk£  +  2s<j>  =  (2q  +  1)tt 


where  q  =  0,  1 ,  2 ,  3,  or,  since  k  =  y-  =  -  ~  , 


(A-8) 


(2q  +  1 )  tt  =  2S(j)  + 


4-rrSf  £ 


(A-9) 


If  we  change  the  driving  frequency  so  that  we  go  through  destructive 


(A-10) 


(A-ll) 


interferences,  then  we  have  for  the  beginning  frequency  f -j : 

4irsf-j  £ 

( 2q  1  +  l)ir  +  2s<J)  +  — - — 


and  for  the  final  frequency 


4  TT  S  f «  £ 

(2q2  +  1 )  tt  =  2s  4>  +  — - — 


Subtracting  (A-10)  from  (A-ll),  we  have 


2(q2  -  q,)i»  -  “r  (f2  -  f,)  • 


Letting  q2  -  q^  =  ^  ancl  ^2  ”  ^1  =  ^  anc*  so^v1n9  f°r  the  wave  velocity 
v,  we  have 


APPENDIX  B 


COMPUTER  PROGRAM  FOR  COMPUTING  VELOCITY 


,GROUP= 


, USER= ■ 


2 

23 


//SONIC  JOB  , 

//  PASSWORO= 

**«JOBPARM  L I NES=5 , CAR DS= 5000, ROOM=B I N4 
***ROUTE  PRINT  RMT26 
//  EXEC  SAS , REG  1 0N=  5 1 2K 
//SYS  IN  OD  * 


, T I ME= ( 5 , 0  ) , CLASS=T , 


NOTE:  SAS  OPTIONS  SPECIFIED  ARE: 
S0RT=4 


OPTIONS  LS=72; 
DATA  SOUND: 
INPUT  Q  F; 
CARDS; 


NOTE: 

NOTE: 


DATA  SET  WORK. SOUND  HAS  108  OBSERVATIONS  AND  2  VARIABLES.  2346  OBS 
THE  DATA  STATEMENT  USED  0.08  SECONDS. 


113 

114 


PROC  PRINT; 

TITLE  OBSERVATIONS  ON  F  AND  Q  FOR  GE( 111); 


NOTE:  THE  PROCEDURE  PRINT  USED  0.17  SECONDS 
AND  PRINTED  PAGES  1  TO  2. 


115 

116 


PROC  MEANS; 

TITLE  STATISTICAL  ANALYSIS  FOR  GE(111); 


NOTE:  THE  PROCEDURE  MEANS  USED  0.13  SECONDS 
AND  PRINTED  PAGE  3. 


117 


PROC  SYSREG; MODEL  F=Q; 


NOTE:  THE  PROCEDURE  SYSREG  USED  0.15  SECONDS 
AND  PRINTED  PAGE  4. 


118 

119 

120 


PROC  PLOT; 

TITLE  PLOT  OF  F  VS.  Q  FOR  GE(111); 
PLOT  F*Q; 


NOTE:  THE  PROCCDURE  ’LOT  USED  0.16  SECONDS 
AND  PRINTED  PAGE  5. 


NOTE:  SAS  INSTITUTE  INC. 

SAS  CIRCLE 
PO  BOX  8000 

CARY,  N.C.  27511-8000 


72 


75 


OBSERVATIONS  ON  F  AND  Q  FOR  CE( 111 ) 


OBS 

Q 

14:14  THURSDAY 
F 

55 

40 

28.144 

56 

50 

28.714 

57 

60 

29.302 

58 

70 

29.865 

59 

80 

30.429 

60 

90 

30.995 

61 

100 

31.564 

62 

110 

32.150 

63 

120 

32.714 

64 

130 

33.296 

65 

140 

33.817 

66 

0 

25.863 

67 

10 

26.425 

68 

20 

26.996 

69 

30 

27.568 

70 

40 

28.155 

71 

50 

28.744 

72 

60 

29.301 

73 

70 

29.865 

74 

80 

30.430 

75 

90 

30.998 

76 

100 

31.581 

77 

1  10 

32. 146 

78 

120 

32.717 

79 

130 

33.303 

80 

140 

33.811 

81 

0 

25.792 

82 

10 

26.373 

83 

20 

26.943 

84 

30 

27.497 

85 

40 

28.068 

86 

50 

28.652 

87 

60 

29.22B 

88 

70 

29.801 

89 

80 

30.376 

90 

90 

30.938 

91 

100 

31.503 

92 

110 

32.085 

93 

120 

32.673 

94 

130 

33.235 

95 

140 

33.801 

96 

0 

25.792 

97 

10 

26.318 

98 

20 

26.889 

99 

30 

27.481 

100 

40 

28.036 

101 

50 

28.610 

102 

60 

29.183 

103 

70 

29.779 

104 

80 

30.318 

105 

90 

30.886 

106 

100 

31.452 

107 

1  10 

32.043 

108 

120 

32.625 

.^>XaV>V 


V'V  j,"1*  n”v  / 


76 


1 


1 

iS 


77 


STATISTICAL  ANALYSIS  FOR  GE ( 111)  4 

14:114  THURSDAY,  DECEMBER  5.  1985 


MODEL:  MODELOI 

SSE 

0.576765 

F  RATIO 

132356.05 

DEP  VAR:  F 

DFE 

MSE 

106 

0.005441183 

PR0B>  F 
R-SQUARE 

0.0001 

0.9992 

VARIABLE 

DF 

PARAMETER 

ESTIMATE 

STANDARD 

ERROR 

T  RATIO 

PROB> I T I 

INTERCEPT 

1 

25.774174 

0.013409 

1922.1649 

0.0001 

Q 

1 

0.056941 

0.0001565134 

363.8077 

0.0001 

78 


PLOT  OF  F  VS.  Q  FOR  GE(111)  5 

14:14  THURSDAY,  DECEMBER  5,  1985 
PLOT  OF  F*Q  LEGEND:  A  =  1  OBS,  B  =  2  OBS,  ETC. 


25  + 

I 

-+- 

0 


- + - + - 

TOO  120  140  160 
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